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Abstract 

We investigate how the pattern of contacts between species in mutuahstic ecosystems is affected by the 
phylogenetic proximity between the species of each guild. We develop several theoretical tools to measure 
that effect and we use them to examine some real mutualistic sytems. We aim at establishing the role of 
such proximity in the emergence of a nested pattern of contacts. We conclude that although phylogenetic 
proximity is compatible with nestedness it can not be claimed to determine it. We find that nestedness 
can instead be attributed to a general rule by which species tend to behave as generalists holding contacts 
with counterparts that already have a large number of contacts. A nested ecosystem generated by this rule, 
shows high phylogenetic diversity. This is to say, the counterparts of species having similar degrees are not 
phylogenetic neighbours. 

Keywords: Nested networks; Mutualistic communities; Phylogenetic proximity; Extinction 



1. Introduction 

A sustainable management of ecosystems as well 
as a proper assessment of the impact of human ac- 
tivity on them can only be achieved with a proper 
understanding of the pattern of the interactions be- 
tween the species. We are here interested in the 
case of mutualistic systems. These usually involve 
groups of animals and plants, helping each other 
to fulfill essential biological functions such as feed- 
ing or reproduction. This is the case of systems in 
which animals feed from fruits while dispersing the 
seeds (seed dispersal networks) or insects feed from 
the nectar of flowers while helping the plant in the 
pollination process {■pollination networks). 

The structure of such systems is described by 
means of an adjacency matrix whose elements rep- 
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resent the absence or presence of an interaction be- 
tween the plant and animal species. This informa- 
tion when concerning mutualistic networks strongly 
indicates that they are not a random collection of 
interacting species, but they instead display a high 
degree of internal organization. A pervading feature 
that has been observed is that the adjacency matrix 
is nested, i.e. if species are ordered by increasing 
number of contacts, those of one species turn out 
to be nearly a proper subset of the contacts of the 
next species in the list [l]. This organization in- 
dicates that the generalists of both type of species 
(i.e. those species that interact with a great number 
of species of the other guild) tend to interact among 
them while there are no contacts among specialists 
(i.e. species that interact with very few species of 
the other guild). 

The nested structure of mutualistic networks has 
been attributed to a number of different causes 
and there is still some degree of controversy about 
the ultimate reasons that make this pattern so fre- 
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quently observed. It is fairly obvious that a detailed 
explanation of the interaction behavior of individ- 
ual species can be of little help to understand such a 
generalized pattern that is found across ecological 
systems of very different sizes and types, and in- 
volving plants of different nature and animals that 
range from insects to birds. 

In Ref.y it is claimed that such order may offer 
some advantage for the robustness of the whole sys- 
tem thus suggesting that systems that are currently 
observed are those that have survived less disturbed 
thanks to its nested structure. Other, more elabo- 
rated theories have been proposed. In Ref.[3| nest- 
edness has been attributed to phenotypic affinity 
between species of different guilds while in Ref.[^] 
an extensive analysis is made concluding that phy- 
logenetic proximity could explain the nested orga- 
nization of contacts of some cases of mutualistic 
systems. A somehow different point of view about 
this work can be found in Ref.|5|. 

It has been customary to consider that the oc- 
currence of some positive statistical correlation is a 
sign of causation for the occurrence of the nested 
pattern of contacts. However, the sole fact that in 
a part of the empirical observations two elements 
appear to be statistically correlated should not be 
taken to mean that one is the cause of the other. 
Such correlation may rather indicate instead that 
both elements are not incompatible, i.e. that they 
do not mutually exclude each other or that they 
stem from a third, common cause. One alterna- 
tive way to search for causal relationships is to ex- 
plore the possible dynamic consequences of some 
assumed interaction mechanism, thus validating or 
falsifying hypotheses concerning possible interac- 
tion mechanisms between the species. In Refs. [6|, 
[7| we have proposed a dynamic model, the Self 
Organising Network Model (SNM), that allows to 
study the contact pattern of a system that is consis- 
tent with some hypothetical interaction mechanism 
between mutualistic species. 

One example of this analysis is given by the 
strong positive correlation found between the 
species' abundance and hence the frequency of in- 
teractions, with the pattern of contacts of some 
species l8| . It has been suggested that locally abun- 
dant species are prone to accumulate interactions 
and conversely rare species are prone to lose them 
[9| , as also suggested by neutral theories [lO| . How- 
ever, a dynamic analysis succeeds in generating a 
realistic distribution of species degree under no as- 
sumptions whatsoever on species abundance. A re- 



verse interpretation may therefore appear also to 
be possible, i.e., that accumulation of interactions 
occurs first, and the resulting higher reproductive 
success leads later to local abundance. With such 
dynamic modeling and by comparison with the ob- 
served natural systems, it is possible to reject or re- 
tain some hypothetical interaction mechanisms as a 
possible cause of the observed order. 

In considering real ecological systems one should 
perhaps expect that a single yet universal cause 
may be the prime responsible of inducing nested- 
ness but that other minor causes may help to shape 
down the observed pattern of interactions of each 
individual case . This state of affairs indicates the 
convenience of finding one or more, model indepen- 
dent tests that can directly be performed on obser- 
vational data to gauge the degree in which an ob- 
served pattern of contacts iJJJ] is governed by some 
presumed cause. 

In the present paper we develop several theoreti- 
cal tools to measure the effect of phylogenetic prox- 
imity in the contact pattern of the system and use 
them to examine some real mutualistic systems. We 
aim in this way at establishing whether such prox- 
imity can be taken to be responsible for the emer- 
gence of a nested pattern of contacts, while some 
other minor causes may explain departures in par- 
ticular circumstances. 

2. Theoretical Background 

Mutualistic systems can be analyzed as bipar- 
tite graphs ^3|. The interaction pattern is usu- 
ally coded into a (rectangular) adjacency matrix in 
which rows and columns are labeled respectively 
by the plant and animal species. Its elements 
Kp^a G {0, 1} represent respectively the absence 
or presence of an interaction (contact) between the 
plant species p and the animal species a. The de- 
grees of each plant or animal species can be ob- 
tained from this matrix as: 



a 



(1) 
(2) 



where Kj is the trasposed of the adjacency ma- 
trix. 

In order to explore the connection between the 
phylogenetic structure of plants and animals and 



the corresponding pattern of contacts we will as- 
sume that a phylogenetic tree has been established 
for animals and plants. All species are therefore 
taken as the leafs or tips of a tree-like diagram that 
represents their evolutionary history. 

A prevailing phylogenetic order can be repre- 
sented as a matrix of phylogenetic distances be- 
tween species. In the literature this is obtained 
from the expected covariance of traits between pairs 
of species through a comparison between real data 
sets and the results of a random walk along the 
tree topology. This way of measuring distances is 
based on the statistical model for the transmission 
of characters along the tree and is therefore model 
dependent. Instead of this we use a more direct, 
unambiguous and model independent way defining 
the phylogenetic distance extracting its value from 
the topology of the phylogenetic tree as the number 
of nodes that have to be explored following the tree 
and starting from either species, until a common 
ancestor of the two species is found. The phyloge- 
netic distance d{k, k') between any two species k, k' 
is then defined as the maximum number of nodes 
that are found starting from k or k'. With this def- 
inition the closest distance between any two species 
is 1. If all species are at a distance 1 then all species 
belong to a so-called star-phylogeny. 

This definition is a simplified version of a more 
elaborate ultrametric [14| distance in which a mean- 
ing is ascribed to the length of the branches of the 
tree by relating it to the time involved in the evolu- 
tionary history. The distance that we have defined 
however satisfies all the requirements of a metric 
function and it also provides values that fully agree 
with an intuitive picture. A small value of d{k, k') 
remains associated to species that share the same 
branching sequence in a common evolutionary his- 
tory and a large value corresponds to species that 
have been separated at earlier stages. 

To test for the effects of phylogenetic proximity 
in the pattern of contacts we will consider static 
and dynamic approaches. Among the former we 
will consider the effects of the ecological similarity 
between species, the robustness of the model and 
the direct correlation between the degrees of pairs 
of species and their phylogenetic distance. Among 
the later we will study the stability of some contact 
patterns when different rules governing the interac- 
tion between mutualistic counterparts are consid- 
ered. 



2.1. Static test I: Correlations with the phylogenetic 
distance 

In order to study the infiuence of phylogenetic 
proximity we investigate the dependence of two pa- 
rameters with the phylogenetic distance. One is 
the ecological similarity Sk,k' and the other is the 
dispersion of the degrees G{k) and G{k') of pairs 
of species fc, k' separated by a given phylogenetic 
distance. 

The ecological similarity is defined 



as 



Sk,k' — 



w, 



MP) 



k.k' 



G{k)+G{k') 



(3) 



In ([3]) the matrix W^, ^, corresponds to the pro- 
jected graphs whose matrix elements measure the 
number of mutualistic counterparts that are shared 
by the species k and k' . The matrix W for plants 
or animals are obtained through W^ = KK^ or 
W^ = K^K. The parameter |3] measures the num- 
ber of mutualistic counterparts shared by the two 
species as compared to the total number of coun- 
terparts of both, in this way, <S'/j_j,/| ,„ax = 1/2. In 
order to consider the dependence with the distance 
we consider the average < S{S) > defined as 



< S{d) > = 



1 



E 



s, 



k,k' 



(4) 



(k,k')\d(k,k')=S 



where Ns is the number of pairs of species of the 
same guild that are separated by a phylogenetic dis- 
tance 6. 

If phylogenetically close species share some com- 
mon feeding or polinization strategy, they should 
also be expected to have similar degrees Ref.[J]. 
The second parameter that we consider is therefore 
the dispersion of the degrees of pairs of species as 
a function of the phylogenetic distance. We define: 



AG^.(^)(^) = 



\ 



^ {k,k')\d(k,k')=S 

(5) 

where Gk is the degree of the species k. The sum 
is extended over all pairs of species (fc, fc') that are 
separated by the phylogenetic distance 6. 

2.2. Static test II: Robustness analysis 

Given the adjacency matrix of a mutualistic sys- 
tem it is possible to get a measure of the robust- 
ness of the system. It has been customary Tg] to 



do this through an attack tolerance curve (ATC). 
The basic idea is that if an "attack" is made by 
which some fraction of the species of one guild is 
ehminated, a number of species of the other will 
become extinct because they are left without their 
mutualistic counterparts. It is of course also pos- 
sible to determine the species that survive to the 
above attack. 

In order to use this concept to study the effect 
of the phylogenetic structure of the species that be- 
long to the mutualist system the ideas involved in 
the ATC have to be generalized. Instead of limiting 
the analysis to the number of species that survive 
or become extinct, one can rather study the phylo- 
genetic relation between the surviving species, by 
keeping track of the average phylogenetic distance 
between them. 

In order to do so we construct a Generalized ATC 
(CATC) in which nestedness and phylogeny are si- 
multaneously considered. Nestedness is taken into 
consideration by orderly eliminating the species of 
one guild according to their degrees while the phy- 
logenetic structure of the partially depleted coun- 
terpart guild is tested through the average phylo- 
genetic distance between the surviving species. 

2.3. Dynamic test: the SNM 

Several arguments have been put forward to sup- 
port that the pattern of interactions between the 
two guilds of a mutualistic network is due to some 
underlying reason. Most of those hypotheses lie 
upon statistically significant, observed correlations 
between some parameters. Correlation is however 
not the same as causation. One way out to eluci- 
date a possible causal link between some hypothesis 
about the interaction mechanisms of the species and 
a given interaction pattern, is to explore eventual 
consequences of such hypothesis through a dynamic 
model. 

The basic idea behind this modeling strategy is 
to verify the consistency of the empirically observed 
contact pattern, with some hypothetical interac- 
tion rule that may favor or hamper the contact 
between mutualist species. We refer to such in- 
teraction mechanism as a contact preference rule 
(CPR). An instability may easily be detected when 
the observed adjacency matrix tends to be dras- 
tically changed if species are allowed to redefine 
their contacts to make a better use of some assumed 
CPR. Within the same reasoning it is also possible 
to check the features of the contact pattern that 
emerges from an initial random adjacency matrix 



if iterated changes of this type are allowed to take 
place. This can give a clue of what one should ex- 
pect to observe in nature if some given CPR is the 
prevailing interaction mechanism among the species 
of system. This kind of analysis allows to trace 
possible causal relationships between a generalized 
individual behavior and the features of the global 
pattern of contacts. 

From a purely theoretical point of view the above 
considerations are equivalent to consider that the 
observed pattern of interactions corresponds to an 
optimal assignment of the contacts between both 
guilds, with two constraints, first the fulfillment of 
some hypothetical, generalized, individual behavior 
trait and second, the total energy expenditure of all 
the species of the ecological system in their search 
for contacts. In other words one may attempt to 
describe the observed pattern of contacts as the re- 
sult of a (combinatorial) optimization problem by 
which contacts in the adjacency matrix are placed 
in such a way as to reach an extreme of some utility 
function defined in terms of the prevailing CPR. 

An exaniple of a dynamic model of this kind is 
the SNM [2] [6|. This has been introduced in or- 
der to account for nestedness in mutualistic webs. 
Within that model, the mutualistic system is as- 
similated to a bipartite graph and the topology of 
the corresponding network is established as the re- 
sult of a self-organization process. This amounts to 
progressively redefine the links obeying some hypo- 
thetical CPR that is assumed toprevail among the 
nodes of the network. In Refs.[4|,[6| we show that 
a CPR where species tend to have contacts with 
counterparts that already hold a greater number of 
contacts leads to nested networks. 

The fulfillment of the additional constraint of 
limiting the total energy consumed by all the 
species belonging to the system is built into the 
model by imposing that the self organization pro- 
cess takes place keeping constant the total number 
of contacts. It is thereforer implicitly assumed that 
on the average, the establishment of each contact 
involves some fixed amount of energy. One could 
consider that this is spent in the process of search- 
ing for the counterpart. 

It is worth to stress that the self-organization pro- 
cess does not represent a real life behavior of plants 
and animals of the system. It therefore does not aim 
at reproducing any evolutive or adaptive process. 
It rather provides a plausible mathematical tool to 
check for instabilities and to search for contact pat- 
terns that are optimal in the sense explained above. 



3. Results: Tests of the effects of phyloge- 
netic proximity 

We now turn to provide a measure of the in- 
fluence of the phylogenetic proximity by focusing 
our attention in a set of systems that have been 
reported in Ref [i] (NCOR, OLAU, ARRl and 
ARR2). Additionally, we take as a reference an 
idealized situation through a benchmark system 
of comparable size that we call PERF (36 plants 
and 50 animals). This is constructed in such a 
way that its adjacency matrix is perfectly nested 
and contacts faithfully follow some assumed (arbi- 
trary) phylogenetic ordering. The phylogenetic tree 
and the adjacency matrix are tailored to perfectly 
match each other, in the sense proposed in [J]. The 
adjaceny matrix together with the (arbitrary) phy- 
logenetic trees are shown in FiglT] 

In all cases we also consider null models consist- 
ing of systems in which the adjacency matrix and 
the structure of phylogenetic trees of both guilds 
are the same as the observed ones but the species 
labeling the tips are randomly permuted. We re- 
fer to these as the " randomized" versions of either 
the real observed system or the ideal PERF-system. 
The null model is constructed considering the aver- 
age over many realizations of such random sorting. 
This is done to provide statistically significant re- 
sults of a situation in which the correlation between 
the structure of the network and the prevailing phy- 
logenetic order has been destroyed. 

3.1. Correlations with the phylogenetic distance 

In Figl2]where we plot the average ecological sim- 
ilarity < S{6) > defined in Eq.((31) as a function of 
the distance i5 separating any pair of species of the 
NCOR and PERF systems. 



It can be seen that for the PERF-system < 
S{6) > decreases monotonically with the phylo- 
genetic distance for plants and animals. This 
monotonous trend is destroyed in the correspond- 
ing randomized system. In the same plots made for 
the NCOR-system there is no qualitative difference 
between the real and the randomized case. 

In FiglSlwe show the plots of AG^(^)(^) as de- 
fined in Eq.® for aU the systems (NCOR,OLAU, 
ARRl, ARR2) together with the one corresponding 
to the PERF benchmark system. If the number of 
contacts were governed by the phylogenetic proxim- 
ity, the smallest value of AG{S) should be found at 
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Figure 1: Adjacency matrix of the benchmark PERF system 
in which a perfectly nested contact pattern is tailored to fol- 
low the (arbitrary) phylogenetic tree of animals and plants. 
Dotted lines are drawn to help the eye in the correspondence 
of contacts with phylogenetic groups 
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Figure 2: Average ecological similarity as a function of the 
phylogenetic distance for plants and animals. We show in the 
same graph the results corresponding to both, the ordered 
and the randomized phylogenetic tree. Upper panel: PERF- 
system. Bottom panel: NCOR matrix. 
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Figure 3: AG as a function of the phylogenetic distance S for 
<5 < 7. Filled (empty) symbols represent animals (plants). In 
every case, the values are normalized by dividing them by the 
average obtained from the null models for the PERF, NCOR, 
OLAU, ARRl and ARR2 systems (AG)j.^^^^^, therefore 
the horizontal full line at 1 represents this value. 



(5=1 because phylogenetically close species must 
have essentially the same num.ber of contacts. Since 
for a greater phylogenetic distance there is also a 
greater diversity of species, AG{6) must also be an 
increasing function of 6. This may no longer hold 
for very large phylogenetic distances because of the 
scarcity of pairs of species in that situation. 

Those two features are clearly displayed in the 
PERF system. Both disappear when the phylo- 
genetic trees are randomized. The curves derived 
from the observed adjacency matrices of the NCOR, 
OLAU, ARRl and ARR2 systems without random- 
ization are compared in Fig|3] to an average of the 
null models for all systems obtained after 100 ran- 
domizations of each phylogenetic tree. They appear 
as random positive or negative fluctuations around 
the randomized average values thus not showing 
any appreciable trend making a difference between 
the observed and the randomized systems. 

3.2. Generalized Attack Tolerance Curves 

We assume that the system is attacked by elimi- 
nating the species of guild A in decreasing order of 
their degrees [l7| . In each stage the GATC that we 
consider provides the average phylogenetic distance 
of the fraction of species of guild B that survive 
|18| . We have found that among a variety of other 
possibilities, this is the plot that more dramatically 
displays the subtle effects of phylogenetic proxim- 
ity. The GATC starts at a value that corresponds 
to the average distance between the species of guild 



B. Since successive stages of the attack produce av- 
erage distances that are either the same or smaller, 
the GATC is supposed to decrease, reaching when 
only one species of guild B is left. 

If contacts were governed by phylogenetic prox- 
imity, close neighbors in the phylogenetic tree of 
guild A must have similar degrees. Assume that 
one starts by eliminating species belonging to one 
group of phylogenetically close generalists of guild 
A. Each elimination stage affects the group of spe- 
cialists of guild B that are in turn assumed to be 
phylogenetic neighbors (see FigH]). Successive ex- 
tinctions of guild A produces little or no change in 
the average distance of survivors of guild B until 
the moment in which the last representative of the 
whole group of specialists of guild B becomes ex- 
tinct. The average distance of the surviving species 
then drops significantly as a whole "branch" disap- 
pears from the tree. This fact leads to "pZateawx" 
in the GATC, putting in evidence the existence of 
groups of phylogenetically close species of guild B 
successively disappearing as the elimination process 
of guild A progresses. 

This GATC can be compared with a null model 
in which the tips of the corresponding phylogenetic 
trees are randomly sorted and many (100) realiza- 
tions of the GATC's are averaged. Upon random- 
izing any trace of phylogenetic proximity is elimi- 
nated but nestedness is nevertheless preserved. The 
GATC is therefore constantly equal to the average 
distance for guild B until the moment in which the 
last species is eliminated and only one species is 
left. 

In FigUwe show several plots of the average dis- 
tance between surviving species as a function of 
the number of elimination steps of their counter- 
parts. In all cases the results of the null model are 
also drawn for comparison. The regular plateaux 
structure can clearly be observed in the benchmark 
PERF system. In the GATC for the NCOR sys- 
tem, traces of phylogenetic affinity can be observed 
both for plants and animals, leading to curves that 
fall well below those of the null model. The GATC 
for the OLAU system only displays some structure 
for plants. In the case of larger systems ARRl or 
ARR2 there is no significant difference between the 
GATC's of the real system and those of the null 
model. 

An interesting consequence of this analysis is that 
the nested structure is not only robust face to at- 
tacks in which concern the preservation of the eco- 
logical system itself but also concerning the preser- 



vation of its ecological diversity. In fact the PERF 
system looses its diversity as measured by the de- 
crease of its average phylogenetic distances well be- 
fore all real systems of comparable size. 



3.3. A dynamic model 

The influence of phylogenetic proximity can be 
cast into the form of a dynamic model of the kind 
explained previously by properly defining a CPR. 
There are two possible algorithms for an SNM ful- 
filling this. They involve the following steps. 

A plant p is chosen at random and elements of the 
adjacency matrix Kp,ao and Kp^ai are also chosen at 



random such that Kp^g = and K,^ 
the following two sums are calculated: 



>-p,ai 



•So = ^d{p,p')Kp^ 



' ,0.0 



p' 



Si = ^d{p,p')Kp,^ai 



1. Next 

(6) 
(7) 



p' 



A swap between this two elements i.e K"^^ — > 1 
and Kp^^ ^> is proposed and it will be accepted 
if, and only if, both animal species ai and oo do not 
become extinct due to the swapping. 

Moreover, one of the two possible extreme CPR's 
defining the model is satisfied: 



• MIN: the condition is 5i > 5*0 

• MAX: the condition is Si < So 

If either of the two conditions is not met, the 
proposed swapping is rejected. In order to proceed 
iteratively, all the above steps have to be repeated 
by interchanging in each step the role of plants and 
animals. 

Within the MIN-CPR possibility, the animal 
counterpart that will finally be selected for p is 
such that all other plants p' that have some con- 
tact with it arc phylogenetically closer to p. This 
is so because the sums 5*0 and 5i involve all phy- 
logenetic distances between the plant p that has 
been selected at random and all other plants that 
make some contact with the two animals oq and oi. 
With this algorithm the configuration of contacts is 
progressively dominated by phylogenetic proximity: 
species of one guild are assumed to interact in the 
same fashion as the rest of species of the same guild 
belonging to their phylogenetic neighborhood. 
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Figure 4: The average phylogenetic distance among surviv- 
ing plants (animals) is plot as a function of the number of re- 
moved plant (animal) species of the other guild. Empty sym- 
bols correspond to the null model (100 random realizations 
of the phylogenetic tree) , and filled symbols to the observed 
systems. Triangles (squares) correspond to plants (animals) 



Within the MAX-CPR possibility, the animal 
counterpart that will be selected for the plant p 
is such that all other plants p' that have some con- 
tact with it are more phylogenetically distant to p. 
In this way the set of species that share contacts 
with p tends to have a greater phylogenetic diver- 
sity. In this alternative all species tend to be as 
generalist as possible as far their phylogenetic clas- 
sifications are concerned. This alternative bears a 
greater similarity with the CPR considered in [6|. 
This is because a larger sum Si or 5*0 is obtained 
not only by involving species that are more phy- 
logenetically distant but also by involving a larger 
number of counterparts. 

In order to check the ordering process generated 
by the SNM we define an effective distance between 
interacting species through: 



D^AP) 



1 j:kM'dik,k')Wk.k' 



(8) 



In Eq.® Wk^k' represents the unweighted adja- 
cency matrix of the projected graphs for animals or 
plants. Its matrix elements are 1 (0) if two species 
share (do not share) mutualistic counterparts. This 
equation provides different results for plants or an- 
imals and should therefore be evaluated separately 
for the two guilds. 

The sum in the denominator of Eq.® is just the 
number of terms appearing in the numerator, there- 
fore D represents the average phylogenetic distance 
between species of the same guild that share at least 
one counterpart of the other guild. D is normalized 
by dividing it by the average phylogenetic distance 
(d) between all plant (animal) species of the sys- 
tem, namely 



{d) = 



Np,a{Np,a - 1) 



(9) 



where Np^A is the number of plant or animal species 
of the system. 

A value D^^^^^ < 1 indicates that phylogenetic 
proximity is a dominant effect, while D^'^"^^ ~ 1 is 
a signature that it is not relevant. We use its values 
to check for the convergence of the ordering process 
implied in the SNM. 

In FigslS] and IH] we plot the values of _D^'(^) 
as a function of the number of iteration steps of 
the SNM for plants and animals for the two possi- 
ble CPR's considered above. Concerning observed 
systems we only consider the NCOR system that 
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Figure 5: Order parameter of Eq.JSjl D^' as a function of 
the number of iteration steps of the SNM algorithm, using 
the MIN-CPR alternative. Filled simbols correspond to us- 
ing as input to the SNM, the observed data of the NCOR 
system while empty symbols correspond to using as input, a 
random adjacency matrix with the same density of contacts 
as the one observed in nature. Triangles (squares) corre- 
spond to plants (animals). 
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Figure 6: Value of the order parameter of Eq. (|8]l D^'^ as 
a function of the number of iteration steps of the SNM al- 
gorithm using the MAX-CPR alternative. Notice that the 
vertical scale is strongly expanded as compared to that of 
FigEl Conventions are the same as in that figure. 



is representative of a possible correlation between 
nestedness and phylogenetic proximity. 

Two different initial conditions are plotted 
in each graph: one corresponds to the adja- 
cency matrix has been reported from observation 
campaigns[J|. The other consists in an adjacency 
matrix of the same number of rows and columns 
and with the same number of contacts except for 
the fact that they are randomly distributed. 

The initial values of D^^^^ are very similar for 
both initial conditions. These are slightly below the 
average phylogenetic distances between species be- 
longing to either of the two guilds. As the number 
of SNM iterations grows both CPR's produce val- 
ues of D^^^^ that reach asymptotic constant values 
that are essentially independent of the initial condi- 
tions, thus indicating that the system has achieved 
a stable and ordered pattern of contacts. However, 
while for the MIN-CPR case it is found that D^^^) 
stabilize at values that are significantly smaller than 
unity, for the MAX-CPR alternative it is found that 
they reach a stable value of 1. 

The initial values oi D^''-^^ indicate that species 
of the same guild sharing at least one mutualistic 
counterpart, are separated by an average distance 
that is very similar to the average distance of both 
phylogenetic trees. The fact that the observed value 
Q^ DP{A) jg QiQsg iQ fji^Q Qjj^g obtained with a random 
matrix suggests that phylogenetic proximity may 
be regarded as an accidental situation involving few 
species of the system. While the MIN-CPR favors a 
contact pattern in which both guilds preserve their 
phylogenetic proximity, the MAX-CPR favors in- 
stead a more diversified situation in which species 
hold contacts with a phylogenetically diverse com- 
munity. 

The SNM algorithm can also be used to test if 
some given pattern of contacts is compatible with 
a CPR involving some kind of phylogenetic depen- 
dence. The asymptotic contact pattern can give a 
clue of what kind of pattern one should expect to 
find for each prevailing CPR. In Figl7]we show the 
adjacency matrix of the SNM system as obtained af- 
ter a great number of iterations of the SNM (panels 
(B) and (C)). These are the asymptotic contact pat- 
terns corresponding to the MIN-CPR alternative 
(panel B) and to the MAX-CPR alternative (panel 
C). In all cases the initial conditions are the em- 
pirically observed contact pattern (shown in Panel 
(A)). These matrices correspond to configurations 
in which D'^^(^) have reached an almost stationary 
value and are therefore nearly optimal in the sense 



explained above. In the same figure the phyloge- 
netic tree of plants and animals are shown to guide 
the eye. 



The MIN-CPR corresponds to a rule in which 
the search of contacts is dominated by phyloge- 
netic proximity. To better understand the emerg- 
ing contact pattern (Panel (B)), one has to bear in 
mind that both animals and plants are considered 
on equal footing. This gives rise to an adjacency 
matrix that breaks in disconnected blocks in which 
phylogenetically close groups of one guild interact 
with similar groups of the other guild. This is the 
opposite of a nested scheme since upon a greater 
number of iterations of the SNM each species tends 
to specialize its contacts as much as possible. By 
the same token, generalists are ruled out of the sys- 
tem. The contact pattern of the NCOR system used 
as an initial condition, becomes therefore severely 
disturbed putting in evidence that it is unstable 
under the presence of the MIN-CPR in which phy- 
logenetic proximity is the dominant rule. 

A similar analysis for the MAX-CPR situation 
shows an opposite behavior. The SNM causes no 
drastic reorderings, reinforcing instead the pres- 
ence of generalists and keeping the matrix mildly 
changed (see panel C). The observed adjacency ma- 
trix must therefore be considered stable under such 
CPR. This run of the SNM also provides additional 
information. The NCOR system hosts a group of 
animals that are phylogenetically close and that are 
all fairly good generalists (the turdus group). Such 
partial correlation between degree and phylogenetic 
proximity is not destroyed by the perturbations in- 
troduced by the SNM, if the prevailing CPR is of 
the MAX type (compare panels (A) and (C)). 

An additional effect of the iteration of the MAX- 
CPR rule is that it leads to an asymptotically sta- 
ble contact pattern that is almost perfectly nested. 
In panel (D) we show the adjacency matrices of 
panels (A) and (C) in which species have been re- 
ordered according to increasing degree, thus com- 
paring the empirically observed nested structure of 
the NCOR system with an asymptotically nearly 
perfect nested pattern produced by the SNM using 
the MAX-CPR. Since the observed NCOR system 
is considerably nested, the effects of the SNM are 
not drastic. 

The occurrence of phylogenetically close species 
with similar degrees and the occurrence of a nested 
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Figure 7; Several adjacency matrices of the NCOR system. Panel(A) the empirical contact pattern with species ordered 
according to the phylogenetic tree (shown along both margins of the matrix); panel(B): contact pattern produced by the SNM 
after 80,000 iterations using as input the empirical matrix shown in panel (A) and the MIN-CPR alternative, species in the 
same order as in panel (A); panel (C): contact pattern porduced by the SNM after 5000 iterations, using MAX-CPR, species in 
the same order as in panel (A); panel (D) same contact pattern as panel (A) (dark pixels) and (C) (hatched pixels) but species 
are ordered by their degree. Dark pixels correspond to observed contacts (Panel (A)) while slanted pattern corresponds to the 
theoretical results (panel (C)) 
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pattern should therefore be considered as indepen- 
dent from each other. When speaking of a cause 
for nestedness one should expect an element that 
is present in all the observed systems, with per- 
haps minor variations in few individual cases. We 
have found here that the only general cause of nest- 
edness is essentially the same as the one reported 
in [6j, i.e. a CPR that tends to place contacts on 
mutualist counterparts that already concentrate a 
large number of contacts. In the present frame- 
work nestedness emerges as a consequence of hold- 
ing contacts with phylogenetically diverse groups of 
species. From this point of view, phylogenetic prox- 
imity has therefore to be considered as compatible 
with a rule inducing nestedenss but is far from be- 
ing a cause of it. 

4. Conclusions 

We have tried several tests to gauge the influ- 
ence of phylogenetic proximity in the occurrence of 
nestedness. They were performed on several real 
systems whose adjacency matrices are reported in 
the literature (NCOR, OLAU, ARRl and ARR2 in 
Ref.[J|) and also on PERF, a benchmark of an ideal 
system in which degrees and phylogenetic proximity 
have been tailored to display a very strong correla- 
tion. 

In these tests we have considered the ecological 
similarity between species, the dispersion in the de- 
grees of species as a function of the phylogenetic dis- 
tances separating them and the shape of the gener- 
alized attack tolerance curves (GATC). In all cases 
we have also compared the results obtained from 
real or ideal systems with those of null models con- 
structed by randomly sorting the tips of the phylo- 
genetic trees. 

The dependence of the ecological similarity or the 
dispersion of the degrees with the phylogenetic dis- 
tance have proven to be essentially undistinguish- 
able from the one obtained for the randomized sys- 
tems in all cases considered. The ideal, perfectly 
organized system shows regular trends of both pa- 
rameters as a function of the phylogenetic distance 
that are not met by real systems. One must how- 
ever bear in mind that both tests are somewhat 
"coarse grained" , in the sense that they involve the 
adjacency matrix as a whole. 

To circumvent this problem we define a GATC 
in which we plot the average distance between the 
fraction of species of one guild that survive after in- 
creasing fractions of the counterpart guild are elim- 



inated. The effect of phylogenetic proximity reveals 
itself as the occurrence of plateaux in an otherwise 
monotonously decreasing function. The results in- 
dicate that the influence of phylogenetic proxim- 
ity is more an accident than a rule. It can only 
be detected in the NCOR system while in the rest 
(OLAU, ARRl and ARR2) only traces can some- 
time be observed. Moreover the results indicate 
that the nested structure of the network makes the 
system robust in terms of ecological diversity. In 
fact, considering the average phylogenetic distance 
as a measure of the ecological diversity, real systems 
preserve this better than the PERF system. 

In order to investigate for possible causes of nest- 
edness rooted in phylogenetic proximity we have 
complemented the static tests mentioned above 
with the use of the SNM. An example of the out- 
come of this test is reported focusing only in the 
NCOR system. 

The general conclusion that stems from the dy- 
namic tests is that an interaction between species 
that exclusively prefers phylogenetic proximity can 
never give rise to a nested contact pattern. A nested 
contact pattern turns out to be unstable in the pres- 
ence of such interaction rule. A way to see this is by 
realizing that such interaction mechanism relies in 
the generalized occurrence of species that are spe- 
cialists therefore ruling out generalists. However, 
these are an indispensable ingredient of a nested 
organization. A contact rule governed by phylo- 
genetic proximity for both mutualists guilds there- 
fore tends to destroy a nested pattern of contacts. 
These tests indicate in addition that the adjacency 
matrices would tend to break down into separate, 
nearly independent components in which groups of 
phylogenetically close neighbors of both guilds hold 
contacts among each other but not with the rest of 
the species of the ecological system. The stabil- 
ity analysis performed using the SNM proves that 
a contact pattern such as the ideal PERF system 
is unstable under the perturbations produced by a 
contact rule favoring phylogenetic proximity. 

We have also proved that an alternative inter- 
action pattern dominated by phylogenetic diversity 
is instead a much better approach to describe real 
situations. This interaction mechanism is one in 
which species are assumed to hold contacts with 
counterparts that are already visited by a greater 
diversity of species. This could be interpreted as a 
stylized version of a rule by which all species tend 
to put the least possible requirements in their feed- 
ing or poUinizing counterparts. It is worth to stress 
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that this rule is fully consistent with the ones tested 
in Refs.[6| and ^ and provides an additional proof 
of the validity of the possible sources of nestedness 
mentioned there. In those situations highly realistic 
degree distribution functions and contact patterns 
are produced by the SNM by only assuming that 
species tend to hold contacts with species that al- 
ready hold a greater number of contacts. 

The second contact rule by which one species 
tends to visit a set of counterparts having a max- 
imal diversity agrees in fact with one of the two 
schemes suggested in the opening figure of Ref . [J| . 
In addition if a group of phylogenetically close 
species happen to have similar contact patterns, the 
corresponding contact pattern turns out to be sta- 
ble under such maximal diversity interaction rule. 
A set of phylogenetically close species of the NCOR 
system that also are good generalists is stable un- 
der the perturbations of the SNM. This points in 
the direction that, although phylogenetic affinity 
can not determine the existence of nestedness it is 
fully compatible with it. 

The results obtained with the dynamic tests are 
consistent with those obtained with the static ones. 
Both point into the same direction of considering 
the correlation of phylogenetic proximity and de- 
gree, if any, is largely accidental and that it can 
never be considered as a general cause of nested- 
ness. A dominant cause of the generalized nested- 
ness found in mutualistic ecosystems perhaps lies 
on the simple fact that species that we observe in 
real systems today are those that tend to put the 
least possible restrictions on their mutualist coun- 
terparts. 



Acknowledgments 

The authors wish to acknowledge helpful discus- 
sions and criticism from D. Medan and M. Devoto. 



References 

[1] Bascoinpte, J. Jordano, P. Melian, C. J. and Olesen 
J.M., 2003. The nested assembly of plant-animal mutu- 
alistic networks. Proc. Nat. Acad. Sci. USA 100, 9383- 
9387; Jordano, P., Bascompte, J. and Olesen, J. M. 
In Plant-pollinator interactions: From specialization to 
generalization N. Waser and J. OUerton, ed.. University 
of Chicago Press. 173-199,(2006) 

[2] Burgos E., Ceva H., Perazzo, R. P. J., Devoto M., 
Medan, D., Zimmermann M., and Delbue A. M., 2007. 
Why nestedness in mutualistic networks, J. of Theor. 
Biol., 249, 307-313 



[3] Jordano P., Bascompte J., Olesen J. M. 2003. Invarinat 
properties in coevolutionary networks of plant-animal 
interaction. Ecology Lett. 6 pp 69-81 
[4] Rezende E.L., Lavabre J.E., Guimaraes Jr. P.R., Jor- 
dano P., Bascompte J., 2007. Non-random coextinc- 
tions in phylogenetically structured mutualistic net- 
works. Nature 448, 925-928 
[5] Susanne S. Renncr 2007. Structure in mutualistic net- 
works. Nature 448, 877-878. 
[6] Medan, D., Perazzo, R. P. J., Devoto, M., Burgos E., 
Zimmermann M, Ceva H. and Delbue A. M., 2007. 
Analysis and assembling of network structure in mu- 
tualistic systems. J. of Theor. Biol., 246, 510-521 
[7] Burgos E., Ceva H., Hernandez L., Perazzo, R. P. J., 
Devoto M., Medan D. 2008. Two classes of bipartite 
networks: nested biological and social systems. Physical 
Review E78, 046113 
[8] Vazquez D.P., Aizen, M.A., 2006. Community-wide 
patterns of specialization in plant-pollinator interac- 
tions revealed by null models. In: Waser, N., OUerton, 
J. (Eds.), "Plant-Pollinator Interactions. From Special- 
ization to Generalization" . University of Chicago Press, 
Chicago, 200-219. 
[9] Stang, M., Klinkhamer, P.G.L., van der Meijden, E., 
2006. Size constraints and flower abundance determine 
the number of interactions in a plant-flower visitor web 
Oikos 112, 111-121. 

[10] D.P. Vazquez, N.P. Chacoff and L. Cagnolo, 2009. Eval- 
uating multiple determinants of the structure of plant- 
animal mutualistic networks. Ecology 90, 2039-2046. 

[11] T. M. Lewinsohn, Paulo I. Prado, P. Jordano, J. Bas- 
comte and J. M. Olesen 2006 Structure in plant-animal 
interaction assemblages. Oikos 113,174-184 

[12] Atmar W. & Patterson B.D., 1993. The measure of or- 
der and disorder in the distribution of species in frag- 
mented habitat. Oecologia 96, 373-382 

[13] Newman M.E.J. , Strogatz S.H. and Watts D.J., 2001. 
Random graphs with arbitrary degree distributions and 
their applications. Phys.Rev. E64, 026128 

[14] R. Rammal, G. Toulouse, M. Virasoro, 1986. Ultra- 
metricity for physicists Rev. Mod.Phys. 58, 765-788 

[15] E. Burgos, H. Ceva, L. Hernandez, R.P.J. Perazzo, 
2009. Understanding and characterizing nestednes in 
mutualistic bipartite networks. Comp.Phys.Comm. 180, 
532-535 

[16] Memmott, J., Waser, N.M., Price, M.V., 2004. Tol- 
erance of pollination networks to species extinctions. 
Proc. R. Soc. B 271, 2605-2611; Memmott, J., Alonso, 
D., Berlow, E.L., Dobson, A., Dunne, J. A., Sole, R., 
Weitz, J., 2005. Biodiversity loss and ecological network 
structure. In: Pascual, M., Dunne, J. A. (Eds.), Ecolog- 
ical Networks: Linking Structure to Dynamics in Food 
Webs. Oxford University Press, New York, 325-347 

[17] This is the only order that allows to see some structure 
in the ATC of nested systems. In the limit of perfect 
nestedness, an attack in increasing order only shows its 
effect when the last species is eliminated. 

[18] As a complement, one can also consider the fraction 
of species of guild B that become extinct {secondary 
extinctions) because they are left without a mutualistic 
counterpart in guild A. 



12 



